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ABSTRACT 

Quillen et al.(2007) presented an imaging survey with the Spitzer Space Tele- 
scope of 62 brightest cluster galaxies with optical line emission located in the 
cores of X-ray luminous clusters. They found that at least half of these sources 
have signs of excess infrared emission. Here we discuss the nature of the IR emis- 
sion and its implications for cool core clusters. The strength of the mid-IR excess 
emission correlates with the luminosity of the optical emission lines. Excluding 
the four systems dominated by an AGN, the excess mid-infrared emission in the 
remaining brightest cluster galaxies is likely related to star formation. The mass 
of molecular gas (estimated from CO observations) is correlated with the IR lu- 
minosity as found for normal star forming galaxies. The gas depletion time scale 
is about 1 Gyr. The physical extent of the infrared excess is consistent with that 
of the optical emission line nebulae. This supports the hypothesis that the star 
formation occurs in molecular gas associated with the emission line nebulae and 
with evidence that the emission line nebulae are mainly powered by ongoing star 
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formation. We find a correlation between mass deposition rates (Mx) estimated 
from the X-ray emission and the star formation rate estimated from the infrared 
luminosity. The star formation rates are 1/10 to 1/100 of the mass deposition 
rates suggesting that the re-heating of the ICM is generally very effective in re- 
ducing the amount of mass cooling from the hot phase but not eliminating it 
completely. 

Subject headings: stars: formation - galaxies: clusters: general - galaxies: active 
- galaxies: elliptical and lenticular, cD - (galaxies:) cooling flows - infrared: 
galaxies 



Introduction 



The hot T~ 10 7-8 K X-ray emitting gas is currently thought to constitute the bulk of the 
baryonic mass in rich clusters of galaxies. An important aspect of the overall physics of the 
intracluster medium (ICM) concerns the central regions of clusters (r < 10 — 100 kpc), where 
the inferred ICM densities and pressures in some cases are sufficiently hig h that cooling to 



T S, 10 4 K can occur on tim e scales shorter th an the cluster lifetime (e.g., ICowie &: Binney 



19771 : iFabian fc Nulsenl 119771 ; lEdge et al.l Il992l ) . These "cooling core" clusters often exhibit 



intense optical emission-line nebulae associated with the centrally dominant (cD) galaxies 
at their centers, together with blue continuum excess emission, and the strength of these 
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The previous paradigm pictured the ICM as a relatively simple place where gas cooled 
and slumped in towards the center of the cluster in a co oling; flow with mass accretion rates 
of hundreds of solar masses per year (e.g., lFabianlll994j ). However, X-ray spectroscopy with 



XMM-Newton and Chandra has failed to find evid ence for gas at temperatures below about 



one-third of the cluster virial temperature (e.g., iKaastra et al.l l200ll ; iTamura et al.l 12001 



Peterson et aPboOll . [iooil : IPeterson & FabianlEooeh . The limits on the luminosity of the in 



termediate temperature gas imply reductions in the inferred mass accretion rates by factors 
of 5-10. Recent theoretical models indicate that intracluster conduction, combined with an 
episodic heat source in the cluster core, such as an AGN or star formation, are candidates for 
explaining both the X-ray emission from cluster cores and the opti cal emission-line phenom- 
ena associated with the cores with these rapid-cooling spectra (e.g., iRuszkowski &: Begelman 



20021 : IVoigt et al.ll2002l : IFabian et al.ll2002 



Narayan &: Medvedev Il200ll ). One widely consid- 



ered possibility is that an important source of heat in the ICM are bubbles driven by radio 
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galaxies (e.g. 


Baum & O'Dea 1991 
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) which halts the cooling of the gas. The ICM now 



appears to be a very dynamic place where heating and cooling processes vie for dominance 
and an uneasy balance is maintained. Since these same processes may operate during the 
process of galaxy formation, the centers of clusters of galaxies provide low redshift laborato- 
ries for studying the critical processes involved in galaxy formation and supermassive black 
hole growth. At the present time, the main questions are (1) How much gas is cooling out of 
the ICM? (2) How much star formation is ongoing? (3) What is the impact of the gas and 
star formation on the central BCG? 

As little mass is needed to power the AGNs at the center of bright cluster galaxies 
the only way to remove cooled gas from the ICM is to form stars. Measurements of the 
star formation rate in cluster galaxies can therefore provide constraints on the efficiency 
of cooling, the fraction of gas that cools and so the needed energy input to prevent the 
remainder of the gas from cooling. It is also possible that the ICM in cluster galaxies is not 
in a steady state or experiences periods of enhanced cooling and star formation and periods 
of relative activity when cooling is prevented. Star formation and associated supernovae 
also provide a source of mechanical energy, tho ugh this is not sufficient to match the X-ray 
radiative energy losses ( McNamara et alTbood ). 



ISO observations detected the cluster Sersic 159-03 ([Hansen et al.ll2000l ). Recent Spitzer 
observations have demonstrated that star formati on is common in cooling core BCGs (jEgami et al 
2006al ; iDonahue et al.l l2007bl ; IQuillen et al.l 120071). An in frared excess is found in about half 
of the sample of 62 BCGs studied b y lQuillen et al.l (120071 ) (Paper I). In this paper we discuss 
the results of IQuillen et al.l (120071 ). We examine correlations in the data and discuss the 
implications for star formation in BCGs and the balance of heating and cooling in the ICM. 
Specifically we search for correlations between star formation rates, radio, Ha, CO and X-ray 
luminosities and mass deposition rates estimated from the X-ray observations. In this paper 
all luminosities have been corrected or computed to be consistent with a Hubble constant 
H = 70 Mpc -1 km s _1 and a concordance cosmology (Qm — 0.3 and flat). 



2. Comparison data 



Quillen et al. 


(2007) 


Quillen et ali 


(2007) 



The properties of the BCG sample are discussed by 
data for the BCGs in our sample are listed in Table 1 by 
able, this Table lists X-ray (primarily ROSAT 0.1-2.4 keV), radio (1.4 GHz), and Ha lumi- 
nosities (from long-slit spectra and SDSS data) and [OIII](5007A)/Hg flux ratios. Brightest 
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cluster galaxies can host both star formation and an active galactic nucleus. X-ray lumi- 
nosities provide a constraint on the mass in and radiative losses from the hot ICM. The Ha 
recombination line is excited by emission from hot stars produced during formation or from 
an AGN. We note emission lines are detected in ~ 10 — 20% of typical optically selected 
BCG s, ~ 30 — 40% of X-ray selected BCGs, and almost 100% for BCGs in cooling core clus- 
ters JPonahue et al.lll992l : [Crawford et aDll999l : Best et al.ll2007i fedwards et aliboOTT ). To 
discriminate between the presence of an AGN and star formation we have sought a measure 
of the hardness of the radiation field through the [OIII](5007A) /H/5 optical line ratio. Fluxes 
in the radio also provide a constraint on the properties of the AGN. Below we discuss star 
formation rates estimated using infrared luminosities derived from aperture photometry also 
listed in Paper I, molecular gas masses estimated from CO observations and mass deposition 
rates measured from X-ray observations. The statistical tests for correlations between the 
various quantities are given in Table [TJ 



Estimated star formation rates 



Bell 


2003; 


Calzetti 


2008) 



to estimate a star formation rate (SFR) (e.g. 
to consider whether some sources h ave a contribution t o the IR from a Type II AGN with 
an optically bright accretion disk. iQuillen et all ((2003) identified Z2089,A1068, A2146 as 
likely to have an AGN contribution based on red 4.5/3.6 micron color, unresolved nucleus 
seen in IRAC color maps and high [OIII]/H/3 flux ratio. R0821+07 was flagged as possibly 
similar as it has an unresolved nucleus in IRAC color map and a high [OIII]/H/3. It also 
has a remarkably red 8.0/5.6 micron color similar to a Seyfert 2 with an embedded dusty 
AGN. In Figure [T] we plot the ratio of 4.5/i and 3.6/zm fluxes against redshift (data from 
Paper I). The clear trend seen is as expected for a passive stellar population but with a few 
notable exceptions. The sources with strong [OIII] (Z2089, A1068, and A2146) lie above the 
trend as do A2055 and A2627, that show evidence for a BLLac continuum in optical spectra 
(jCrawford et al.lll999l ). The two galaxies that lie below the trend are Z2072 and Z9077 and 
are among the fainter objects in our sample. Z9077 was the only object not detected at 24 
microns. The four that are between z=0.09 and and z=0.15 and lie slightly above the trend 
are A1885, A2055, A2627 and R0352+19. It's not obvious why these 4 sources lie above the 
trend or if this is significant R0352+19 and R0821+07 are quite red in the 8/5.8 micron color, 
and R1532+30 and Z348 are pretty red in 8/5.8 but do not stand out in the 4.5/3.6 color vs. 
z plot. Thus, the combination of diagnostics (4.5/3.6 micron color, red unresolved nuclear 
source, and high [OIII] /H/5 ratio) identifies some sources with a strong AGN contribution. 
The remaining objects are likely to be free of strong AGN contamination. 
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Previous optical and UV observations have found evidence for s ignificant star for 
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Raffertv et all2006h . Table 1 of buillen et all (120071 ) lists [OIII]/H/3 



199 



Cardie 



lO'Dea et al. 



et al. 



2004 



ratios for most of the BCGs. Except for the few which may host an AGN, the ratios are 
consistent with the gas being ionized by hot stars. 



In Table [2] we present the estimated infrared luminosities from iQuillen et al.l (120071 ) 
and the estimated star formation rates. The SFR rate can be estimated from the infrared 



luminosity with equation 5 by iBelll (120031 ): 



V(M yr" 



.4 



L 



IR 



L 







IRJ 



Here the constant A = 1.57 x 10 -10 for L IR > 10 n L o and A = 1.17 x 10 -10 at lower 
luminosities. The SFRs are in the range of about 1 to a few tens of M yr _1 . The objects 
with SFR above about 50 M yr" 1 are likely AGN dominated. 

In Table [3] we list available SFRs in different wavebands. We see that there is dispersion 
in the estimated SFRs. However because of the effects of dust and geometry we do not 
necessarily expect agreement between SFRs estimated in the IR vs. the UV/optical. Much 
of the variation can be accounted for by aperture mismatch, differences in assumptions 
about star formation history i.e., burst vs. constant star formation, extinction, and perhaps 
differences in the amount of dust available to re-radiate in the FIR. Note that Abell 1068 and 
A2146 show large discrepancies between our FIR SFR and the U-band SFR and both are 
flagged as possible AGN. Given the expected dispersion, the rough agreement between the 
star formation rates is consistent with the IR emission being dominated by star formation. 



3.1. Caveat - the Dust-to-Gas Ratio 



The derived SFR might be underestimated if the cold gas in the BCGs has a low dust- 
to-gas ratio. This might be the case if the gas has cooled from the hot ICM and if the dust 
was destroyed while in the hot phase and there has not been sufficient time to form dust at 
the levels typically seen in normal star forming galaxies. However, t here are several argu- 
ments against a low gas-to-dust rati o . (1) The observations of H 2 (e - g., 



Edge et al.l [2OO2I: lHatch et al 



2007J) and CO (lEdge 



2001 



2005 



Jaffe et al 



Salome fc Combesi 



2005 



Donahue et al 



Egami et all l2006bl : 



2000: 



Johnstone et al. 



20031 . 120041 ) associated with the BCG opti- 



cal emission line nebulae require the presence of significant amounts of dust to shield the 
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ters (e.g., 


Sparks et al. 


Koekemoer et al. 


1999). 



1999h . (3) Studies of the nebulae in the BCGs of cool core clusters suggest 



the presence of dust-to -gas ratios consistent with Galactic values (jSparks et al.lll989l . Il993 



Donahue fc Voitl I1993T). (4) Theoret i cal arguments suggest that dust could form quickly 



inside cool clouds (IFabian et al.lll994j ; IVoit fc Donahueill995l ). 



4. Comparison between infrared luminosity and X-ray luminosity 

The integrated X-ray luminosity of a cluster is dependent on the combination of its core 
and larger scale structure. As such any correlation between this global property and the 
properties of the BCG may indicate an underlying link, particularly as our sample from its 
selection will favor cool cores. Therefore we plot the X-ray luminosity of the host cluster 
(listed in Table 1 of Paper I) against estimated infrared luminosities for all BCGs with color 
ratio Fgfjja/ F 5 Sflia > 0.75 in Figure [2J In Figure [3] we show X-ray luminosities compared 
to the color Js /tm / J fj.s/i.nr This study covers a much larger range in X-ray luminosity than 



Egami et al.l (l2006al ). We see that BCGs with higher IR luminosity and redder 8 to 5.8/im 
colors (indicating an IR excess) tend to have higher X-ray luminosities. Though, there are 
many objects with high X-ray luminosity (Lx > 10 44 erg/s), which do not have an IR excess. 

It is interesting to compare the kinetic energy injected by supernovae (from a star 
formation rate consis tent with the infrared luminosity) to the energy radiated in X-rays. 



Leitherer et al.l (119991 ) estimate a mechanical energy of about 10 42 erg s _1 normalized for a 
star formation rate of 1M Q yr _1 . These conversion factors have been used to estimate the 
mechanical energy due to supernovae as a function of infrared luminosity. This relation is 
shown as a dashed line in Figure [2j We see that there are a few BCGs for which there may 
be sufficient mechanical energy to resupply the X-ray luminosity. Howev er, in general, for 



the s ample as a whole, we confirm the finding of previous studies (e.g., iMcNamara et al. 



20061 ) that mechanical energy input from supernovae is not sufficient (by a few orders of 
magnitude) to account for the current radiative energy losses of the intracluster medium as 
a whole or in the core. 



5. Comparison to radio luminosity 



We find a modest (almost 3a) correlation between the infrared luminosity and the radio 
luminosity at 1.4 GHz (as we show in Figure H]). We compare the radio fluxes to those 
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appropriate for star forming objects with a solid line on the lo wer right in Figure HI The 
radio-IR relation for star forming objects (equation 3; iBelll 120031 ) 



erg cm -2 s _1 Hz _1 



3.75 x 10 12 +9 erg cm^ 2 s" 1 Hz" 1 



(2) 



where q is a logarithmic index. We have used the mean value q = 2.34 by lYun et al.l (]200ll ). 



On Figure HI the majority of radio fluxes are well above this relation. The three objects 
below the line are (from left to right) NGC4104, R0821+07, A1068. NGC4104 is nearer 
than the other objects in the survey and it is possible that the Ha flux and radio flux 
density have been underestimated. The other two clusters (Abell 1068, R0821+07) have 
F8fj, m / F 5 S ^ m > 1.3 and have unresolved red sources seen in the IRAC color maps and so are 
likely to be dominated by an AGN. 

Thus, the BCGs (independent of whether they have an IR excess) tend to have radio 
emission which is dominated by that produced by an AGN. Based on hot IR colors and high 
[OTII]/H/3 ratio s it appears that o nly 4 of the BCGs host a Type II AGN with a luminous 
accretion disk (jQuillen et al.ll2007l ). Thus, either the AGN in most of the BCGs are currently 
turned off, or they are accreting in a low luminosity mode. The (weak) correlation between 
radio and IR luminosity may be a consequence of the correlation between mass accretion rate 
and SFR §[5} i.e., the cooling gas feeds the AGN and makes gas available for star formation. 
In addition, the ratio of mechanical energy in the radio source o utflow to the radio luminosity 
can vary by about 3 orders of magnitude iBirzan et al.l (120041 ) . Thus, the radio luminosity 
can be a poor measure of the impact of the radio source on its environment. 



6. Comparison to Ha luminosity 



We compare the Ha luminosities from limited aperture spectroscopy to the infrared 
luminosities in Figure [5] finding a strong correlation between the two. We also see a cor- 
relation in Ha flux vs. 24 micron flux density (Figure [6]). These correlations show that 
the Ha and infrared emission arises from the same or a related power source. We suggest 
that the dominant power source for the Ha and infrared emission is star formation. This is 
consistent with previous evidence that the optical emission line nebulae are mostly powered 



by UV phot ons from young stars with a possib 



amsm (e.g., 



Johnstone fc Fabian 



1999; O'Dea et al. 



2004 : 



1988 



Wilman et al 



Allen 



2006 



e sec o ndary contribution from another mech- 



1995 



Hatch et al 



Voit fc Donahue 1997; Crawford et al. 



20071 1 The Ha-SFR law relating 



the star formation rate to the Ha luminosity 

SFR (M yr- 1 ) = 



L(Ha) 



1.26 x 10 41 erg s" 1 



(3) 
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( IKennicuttl Il998l ) is shown as a dashed line on the plot. We have scaled the line down 
by a factor of 2.8 because our Ha luminosities are uncorrected for reddening. We see in 
Figure that the points tend to lie a factor of a few below this line, i.e., the observed lu- 
minosity in Ha is lower than that expected from the estimated infrared luminosity. The 
discrepancy is larger at lower X-ray luminosity Lx < 10 44 ergs/s. Our Ha luminosities 
are taken mainly from spectroscopy with a long slit of width l'.'S (jCrawford et al.l Il999l ) 
or the 3" diameter fibers of the SDSS. Narrow band Ha+[NII] images and IFU observa- 
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: I 


iatch et al. 
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). We find that the nebu- 



lae are all larger than the spectroscopic apertures, with a median size of 7"1 (geometric mean 
of major and minor axes). The total fluxes determined from the narrow band images and 
IFU spectroscopy are larger than those from lCrawford et al.l (119991 ) or SDSS, with a median 
ratio of 1.4. Thus, it seems likely that aperture effects contribute to the Ha deficit, though 
larger samples with narrow band imaging are needed to determine whether this can explain 
the whole effect. Additionally, strong absorption of the Ha (relative to normal star forming 
galaxies) could also contribute to the Ha deficit. However, the possibility remains that star 
formation is not the only power source for the Ha and IR emission. If this turns out to be 
the case, it would suggest that there is an additional source of energy which heats the dust, 
but does not i onize the gas. Such an energy source would help to explain the ob s erved opti- 



cal line rat ios (jVoit fc Donahue! 1 199 71 ) and bright H 2 emission (jEdge et al.ll2002l ; iJaffe et al. 
200lLl2005h . 



7. Comparison to molecular gas mass 



We have compiled molecular mass data from lEdgd (120011 ). ISalome &; Combed (120031 ). 
and Edge, in preparation. This subsample consists only of objects that have been surveyed 
for and detected in CO (1-0). The inferred molecular gas masses range from ~ 10 9 to 
~ 10 11 M . Spitzer IRS spectra of the star forming BCG in Z3146 de tect strong molecula r 
hydrogen lines from warm H 2 with an estimated mass of ~ 10 10 M ( lEgami et al.ll2006bl ). 
This provides confirmation that the molecular gas masses can be very large in these BCGs. 

We note that a correlat ion be t ween integ rated molecular gas mass and Ha luminosity 
in BCGs has been found by lEdgd (120011 ) and ISalome fc Combed (120031 ). In Figures [7] and 
[HI we plot the molecular mass against our estimated infrared luminosity and star formation 
rate (listed in T able El) . As found in normal star forming galaxies (e.g., lYoung et al.l 11986 



Kennicuttl 119981 ) . we see a correlation between measured molecular gas mass and both the 
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IR luminosity and the star formation rate in the BCGs. The ratio of molecular gas mass to 
SFR gives a gas depletion time scale which is roughly 1 Gyr. The gas depletion time scale is 
roughly constant over a range of two orders of magnitude in molecular gas mass and SFR. 
Our value of ~ 1 Gyr is i n good agreement w ith the mean value of ~ 2 Gyr found in normal 
star forming galaxies by lYoung et al.l (119861 ) which have molecular gas masses in the range 
(~ 10 9 to 10 10 M ). The long life time of the molecular gas in these BCGs is in contrast to 
the much shorter cooling times for the gas over a range of t emperatures. The hotter phases 
cool on times of ~ 10 6 — 10 8 vr s (IPeterson fc Fabian! 120061 ). while the molecular gas cools 



on even shorter time scales (e.g.. |Jaffe et al.ll200ll ). Given that clusters are relatively young 



(perhaps 4-6 Gyr since last major merger) it is possible that there may have been insufficient 
time for a complete steady-state (cooling leads to cold gas leads to star formation) to be set 
up. 



7.1. The Size Scale of the Star Formation and its Relation to the Optical 

Emission Line Nebulae 

In nearby galaxies there is an empirical relation between star formation rate per unit 
area and molec ular gas surface d ensity. This relation can be described in terms of a Schmidt- 



Kennicutt law (IKennicuttl 1 1 9981 ) 



^ 2) = 2-5 x 10- 4 ( T ^^) ^ , (4) 

where S gas is the surface density of molecular and atomic gas and Esfr is the star formation 
rate per unit area. We can use this relation to estimate the size scale of the star forming 
region. We make the assumption that the star formation is distributed in a region of area 
d\ pc where <ik pc is a diameter in kpc and the surface density S gas = Mn 2 /d%. pc and where M H2 
is the molecular gas mass. Applying this to equation HJ we find a relation between molecular 
mass and star formation rate; 

We have shown this line in Figure [8] computed for diameters d^ pc = 5, 15, and 50. We see that 
the data are consistent with a Schmidt-Kennicutt law, but the diameter of the star forming 
region is not well constrained. The diameter of the star forming region tends to be larger 
for more luminous objects as expected if the diameter is proportional to Lx/Mx- Previous 
studies have shown that the Schmidt law predicts the star formation rates within a factor 
of a few for galaxies over a wide range of morphologies and star formation rates, including 
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starbursts galaxies (IKennicuttl Il998l ) . The previous study of 2 BCGs by (IMcNamara et al. 
20061 ) suggested that the star formation law holds even in BCGs. However, Figure [7| shows 



that at high molecular gas masses Mh 2 > 10 10 M , some BCGs show inferred diameters 
J> 50 kpc which are much larger than suggested by the sizes of the emission line nebulae. 

The sources with the largest estimated star formation regions are R1532+30, A1664, and 
Z8197 with estimated star formation region size scales of 70, 50 and 30 kpc, respectively, 
estimated using the Schmidt type star formation law. None of these is well resolved, all 
have FWHM near the diffraction limit of 7" at 24 microns. For R1532+30 at z = 0.36, the 
FWHM corresponds to a size of 35 kpc. This is below the estimated size of the star forming 
region, R ~ 70kpc. Likewise for A1664, and Z8197, with redshifts of 0.128, and 0.114, the 
FWHM corresponds to about 15 kpc and again this exceeds the estimated size scale by a 
factor of 2-3. The size scale estimates using the star formation law are probably a factor of 
a 2-3 too large for these objects. Those with the smallest estimated star formation regions 
are A85, A262, A2052, and NGC4325 with estimated regions of smaller than 5 kpc. At a 
redshift of 0.0551,0.0166, 0.0351, 0.0259, 7" (diffraction limit at 24 microns) corresponds to 
7, 2.3, 5 and 3.6 kpcs, respectively. The objects with the smallest estimated regions are 
the nearest and so can be resolved in the IRAC images. For A85, the star forming region 
could be the unresolved source at 8 microns that is south east of the brightest cluster galaxy 
nucleus. The brightest cluster galaxy is resolved at 8 microns. For A262 and NGC4235, the 
brightest cluster galaxy is the source of the 24 micron emission and is resolved both at 24 
and 8 microns, consistent with the estimate for the star forming region size of a few kpc. For 
A2052 the brightest cluster galaxy also hosts star formation in its nucleus. The emission is 
unresolved at 24 microns but resolved at 8 microns. This is consistent with the estimated 
size of the star forming region of a few kpc. Except for the case of A85 the estimated sizes 
of the star forming regions of a few kpc are consistent with the sizes estimated from the 
images. In summary, the Kennicutt-Schmidt law gives sizes which are generally consistent 
with those estimated from the images for the small and average sizes, though the largest 
sizes seem to be too large by factors of 2-3. 

We note that the emission line ne bula in cool cores tend to have a bright central region 
with a diameter of order 10 kpc (e.g ., iHeckman et al.lll989l ). with fainter gas extending to 



larger scales (e.g., iJaffe et al.l 120051 ) which is comparable to the inferred size of the star 
formation region. Observations of extended HI absorption in the emission line nebula of 
A2597 suggests that the optical n ebulae are photon bounded and are the ionized skins of 



been found associated with emission line filaments in some BCGs (e.t 



Edge et al.l 12002 ; lHatch et al.l 120051 : IJaffe et al.l 120051 : lEgami et al 



g-, 


Donahue et al. 


200d; 


2006bl: 


Johnstone et al. 



20071 ). Interferometric CO observations show molecular gas associated with the emission 
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line nebula in A1795 (e.g., ISalome fc Combed 12004 ). Also, HST FUV images show FUV 
continuum from yo ung stars associated with the emission line nebulae in A1795 and A2597 
( jO'Dea et al.l 12004 ) . Thus, the spatial association of the FUV, the CO, the H 2 , and the 
optical emission line nebulae suggest that the star formation occurs in molecular gas which 
lies in the optical emission line nebulae. 

Our estimated size scale of d ~ 15 kpc for the star formation region could be biased. 
Objects that have larger and more diffuse star formation regions would have had larger 
molecular gas masses and so would have been detected. Similarly, BCGs with lower and 
more concentrated star formation regions might have been missed. 



8. The Connection between Star Formation and the Properties of the Hot 

ICM 

We use archival Chandra and XMM-Newton observations to calculate X-ray inferred 
mass deposition rates and cooling times for 14 of the selected clusters. We required at least 
15,000 counts from the source in each observation to generate reliable deprojected spectra. 
This restricted the cluster sample to 11 with suitable Chandra archive data and 3 with 
XMM-Newton archive data (R0338+09, R2129+00 and Abell 115). 

The Chandra data were analyzed using CIAOv4.0 beta 2 with CALDBv3.4.1 provided 
by the Chandra X-ray Center (CXC). The level 1 event files were reprocessed to apply the 
latest gain and charge transfer inefficiency correction and filtered for bad grades. Where 
available, the improved background screening provided by VFAINT mode was applied. The 
background light curves of the resulting level 2 event files were then filtered for periods 
affected by flares. For the nearer clusters (z < 0.3), background spectra were extracted from 
blank-sky background data sets available from the CXC and cleaned in the same way as the 
source observations. The normalizations of these cleaned background files were scaled to 
the count rate of the source observations in the 9-12 keV band. For more distant clusters, 
background spectra were extracted from suitable, source-free regions of the source data sets. 

The XMM-Newton MOS data were reprocessed using the emchain task from XMM- 
Newton SAS v7.1.0 to generate calibrated event files from the raw data. Cosmic ray filtering 
was applied by selecting only PATTERNS 0-12 and bright pixels and hot columns were 
removed by setting FLAG==0. Soft proton flares were removed by generating a light curve 
for photons of energy > lOkeV, where the emission is dominated by the particle-induced 
background, and rejecting high background periods. Blank-sky background spectra were 
produced using the blank-sky background event files available from the XMM-Newton Science 
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Operations Center and calibrated according to iRead Ponman I (120031 ). 



Spectra were extracted in concentric annuli centered on the X-ray surface brightness 
peak with a minimum of 3000 counts in each annulus. Chandra spectra were analyzed in 
the energy range 0.5-7.0 keV and XMM-Newton spectra in the range 0.3-10 keV. Suitable 
response files (ARFs and RMFs) were calculated and grouped together with the source 
spectrum, binned with a minimum of 30 counts. 

Deprojected tem perature and density profiles were calculated using a Direct Spectral 
Deprojection method (jSanders fc Fabian 1120071 : ISanders et al.llin preparation! ), which creates 
'deprojected spectra' using a model independent approach, assuming only spherical geome- 
try. Instead of correcting for projection by combining a series of models, this new method 
subtracts the projected spectra from each successive annulus to produce a set of deprojected 
spectra. 



The resulting deprojected spectra were analyzed in XSPEC 11.3.2 (lArnaud Ill996l ). Gas 
temperatures and densities were found by fitting each spectrum with an absorbed single- 
temperature plasma model (phabs(mekal)). The redshift was fixed to the value given in 
Table 1 of Paper I and the absorbing column density was fixed to the Galactic values given by 
Kalberla et al.l (120051 ). The temperature, abundance and model normalization were allowed 
to vary. We used the deprojected temperature and density to determine the cooling time of 
the gas at each radius. The cooling radius was defined to be the radius within which the gas 
has a cooling time less than 7.7 x 10 9 yr, corresponding to the time since z = 1. 

In Figure [9] we plot the Infrared Luminosity vs. the cooling time at a radius of 30 kpc. 
We see th at BCGs with shorte r cooling times have higher IR luminosities consistent with the 
results of lEgami et al.l (j2006al ). This is consistent with the hypothesis that the clusters with 
shorter cooling times have higher star formation rates which result in higher IR luminosity. 

We calculate two different measures of the mass deposition rate. Mi is a maximum 
mass deposition rate, calculated from 



Mr 



2fim H L(< r cool ) 
5k B T(r cool ) 



(6) 



where [imu is the mean gas mass per particle and the luminosity was determined over the 
energy range 0.01-50 keV. L(< r coo i) is directly proportional to the energy required to offset 
the cooling and Mj is a measure of the mass deposition rate if heating is absent. 

Ms was calculated by repeating the spectral fitting to the annuli within the cooling 
radius but now adding a cooling flow model (mkcf low) to the absorbed single-temperature 
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model (phabs (mekal+mkcf low)). The XSPEC model mkcflow models gas cooling between 
two temperatures and gives the normalization as a mass deposition rate, Ms- For each 
spectrum, the temperature of the mekal component was tied to the high temperature of 
the mkcflow component and the abundances of the two components were tied together. 
The low temperature of the cooling flow model was fixed to 0.1 keV. Ms is a measure of 
the maximum rate (upper limit) that gas can be cooling below X-ray tempe ratures and be 



consi s tent with the X-ray spectra. Detailed spectra of nearby bright clusters (IPeterson et al. 



20031 ; iFabian et al.ll2006l ) tend to show an absence of the X-ray coolest gas and indeed for 



the inferred M to be a function of temperature within a cluster. Better quality data for 
the objects here may lead in some cases to lower estimates of Mg- We have listed the mass 
deposition rates for each cluster in Table H] and plotted these values against star formation 
rates estimated from the infrared luminosity in Figure (TUl 



8.1. Implications 



We see that the SFR is proportional to (but significantly less than) the two estimates 
of mass accretion rate. The results show that the star formation rate is about 30-100 times 
smaller than Mj, and 3-10 times smaller than Ms- The observed trends between cooling time 
and the IR luminosity and between Ms and the infrared star formation rates are consistent 
with the hypothesis that the cooling ICM is the sour ce of the gas which i s forming stars. 
Using a nearly identical approach to the X-ray analysis, iRafferty et al.l (120061 ) found a similar 
trend using optical- UV star formation rates. Star formation in these systems is generally not 
heavily obscured, and the star formation rates are approaching and in some cases agree with 
the t he cooling upper lim i ts, Ms- This is cons i stent with the results from X-ray spectroscopy 



Kaastra et al. 2001; Tamura et al. 2001; Peterson et al. 2001. 2003; Peterson fe Fabian 



20061 ) which suggests that most of the gas with short cooling time at a few keV does not 
cool further. Sensitive, high resolution X-ray spectroscopy should soon detect the cooling at 



the level of sta r formation in the Fe XVII lines if this picture is correct (e.g., iSanders et al. 
in preparation! ) . Our fraction of gas which does cool is a mean number and could be affected 
by our Ha selection but we believe the use of a complete X-ray sample will allow this effect 
to be quantified. Nevertheless this number could provide a constraint on the efficiency of 
feedback models that prevent the bulk of the ICM from cooling. If star formation is the 
ultimate sink for the cooling gas, then the fraction of the few keV gas which does cool all the 
way down should be comparable to the ratio SFR/Mx - which we find to be roughly a few 
percent. This fraction is comparable for all the clusters. This suggests that the re- heating 
mechanism (whatever it is) is very effective over a range of size scales and operates nearly all 



the time (i.e., with a short duty cycle) (see iMcNamara fe Nulsen 1120071 ; IPeterson fe Fabian 
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20061 for reviews) . 



8.2. Alternative Energy Sources 



We have proceeded with the assumption that the infrared emission is solely due to star 
formation. Here we examine whether there are reasonable alternative sources of energy for 
this emission. First we consider the hot gas, sin ce dust mixing with the gas can be heated 
and become a source of mid infrared radiation (IDwek Ill986l ; iDwek et al. Ill990l ). We then 
consider cosmic ray and other heat sources. 

The hot gas is potentially a rich energy source which could heat the dust. A consequence 
of such heating is the the energy loss from the gas which means the gas will cool, perhaps 
even exacerbating the cooling flow problem. It offers a solution to the problem seen in the 
X-ray spectra of coo l core clusters in wh i ch gas is observed t o cool down to only about 
1 ke V but no lower (IPeterson et al.l l200ll ; iTamura et al.l l200ll ; iPeterson et al.l 120031 ) . The 
temperature profiles in clusters mean that the coolest gas is at the smallest radii so if there 
is non-radiative cooling of gas at those radii, say due to mixing with cold dusty gas, X-ray 
spectra of the whole core would imply a cooling flow going down to just 1 keV and appearing 
to stop, more or less as obs erved. This can be seen as 'the missing soft X-ray luminosity' 
problem (IFabian et al.ll2002l ) . What is meant by missing soft X-ray luminosity is the emission 
missing from a complete cooling flow when it appears from X-ray spectra to stop at say 1 
keV. 

Figure [TT] shows the missing soft X-ray luminosity for our objects plotted against the 
IR luminosity as open circles. This was obtained by fitting the spectra with a cooling flow 
model which has a lower temperature limit of 1 keV. The missing soft X-ray luminosity is 
then the rate of energy release as that gas cools further to zero K in some non-radiative 
manner. There is a correlation, but the normalization misses by about a factor of 5. This 
means that on average there's 5 times more than needed to account for non-radiative 
cooling of the gas below 1 keV. 

Better agreement can be obtained by increasing the lower fitted temperature above 
1 keV, but in that case the mass cooling rates rise from the more modest rates comparable 
to Ms in Fig. [TD]to the higher, pre-XMM/Chandra values of Mj. This is just because L IR 
is similar to L coo \, the luminosity of the cooling region in the core (i.e., where the radiative 
cooling time is less than say 5 Gyr), which is shown by the filled circles in Fig. [TTJ 



The result is that if dust mixing in hot gas is the only source of infrared emission then we 
have to face cooling rates much higher than can be accommodated in terms of the observed 
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molecular gas ( Edge! I2OO1I ) or observed star formation rates (Table [3]). More plausibly, hot 



gas mixing with dusty cold gas is the source of 10-20 per cent of the infrared emission. In 
this case our results allow for modest rates mass cooling rates of up to tens to hundreds 
M Q yr _1 comparable to the range shown in Table [31 

Cosmic rays also fail as an energy source, unless they are recycled. Since Lir ~ L cool (to 
within about a factor of 3, see Fig fTTl) . then the energy required for the infrared is comparable 
with the thermal energy within r coo i. Consequently the cosmic ray pressure would need to 
be high with a pressure P CR = /cR-Prh with fcR > 0.3 and thermal pressure P-ph- This is 
contrary to the quasi-hydrostatic appearance of the intracluster medium in cluster cores. 

Only if there is some efficient mechanism for energy to flow from the central accretion 
flow / AGN to the dust can an alternative be viable. In the absence of any such mechanism, 
we conclude that the UV radiation from massive star formation must be the energy source 
for the mid-infrared emission measured by Spitzer. 



9. Summary 



Quillen et al.1 (120071 ) obtained Spitzer photometry of a sample of 62 BCGs in X-ray 
bright clusters selected on the basis of BCG Ha flux which tends to favor cool core clusters. 
They showed that at least half of the BCGs exhibit an infrared excess with a luminosity 
Z-ir ~ 10 43 — few xlO 45 ergs/s. In this paper we examined correlations in the data and 
discussed implications for cool core clusters. 

BCGs with an IR excess are found mainly in clusters at high X-ray luminosity (Lx > 10 44 
ergs s^ 1 ). But not all high Lx clusters have a BCG with an IR excess. 

The IR luminosity is proportional to the Ha luminosity, suggesting that they are pow- 
ered by the same or a related source of energy. We suggest that star formation is the 
dominant power source for the the IR and Ha emission. The Ha luminosity falls below the 



Kennicuttl (119981 ) relation probably because the spectroscopic apertures exclude much of the 
extended emission line nebulae. The inferred star formation rates estimated from the IR 
luminosity are in the range about 1 to 50 M yr _1 . In most BCGs, supernovae produced 
by star formation with this SFR cannot account for the X-ray luminosity and so cannot be 
responsible for re-heating the ICM. 

The radio emission in the BCG is dominated by that produced by an AGN rather than 
star formation. However, there is a modest correlation between radio and IR emission. This 
suggests the feeding of the AGN and the fueling of the star formation may have a common 



-17- 



origin, perhaps gas cooling from the hot ICM. 

The mass of molecular gas (estimated from CO observations) is correlated with the IR 
luminosity as found for normal star forming galaxies. The gas depletion time scale is about 
1 Gyr. Given that clusters are relatively young (perhaps 4-6 Gyr since last major merger) 
then it is possible that there may have been insufficient time for a complete steady-state 
(cooling leads to cold gas leads to star formation) to be set up. 

We fit a Schmitt-Kennicutt relation to the molecular gas mass vs. SFR and estimate 
a rough star forming region diameter. For most BCGs the implied sizes of 10-20 kpc is 
comparable to that of the color variations seen in the IRAC data and to the optical emis- 
sion line nebulae. This is consistent with the hypothesis that the star formation occurs in 
molecular gas associated with the emission line nebulae and with evidence that the emission 
line nebulae are mainly powered by ongoing star formation. 

BCGs in clusters with shorter cooling times for the hot ICM have higher IR luminosities. 
We find a strong correlation between mass deposition rates (Mj) estimated from the X-ray 
emission and the SFR. The star formation rate is about 30-100 times smaller than Mj - the 
mass accretion rate derived from imaging, and 3-10 times smaller than Ms - the rate derived 
from spectroscopy. The observed trends between cooling time and the IR luminosity and 
between M$ and the infrared star formation rates are consistent with the hypothesis that 
the cooling ICM is the source of the gas which is forming stars. The correlation between 
mass deposition rates estimated from the X-ray radiative losses and the star formation rates 
suggest that the fraction of gas that does cool is set by the balance of heating and cooling 
by the cooling flow. The low value of SFR/(Mx) suggests that heating is likely to be very 
efficient in preventing most of the gas a temperatures of a few keV from cooling further. 

This work is based in part on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology under 
a contract with NASA. Support for this work at University of Rochester and Rochester 
Institute of Technology was provided by NASA through an award issued by JPL/Caltech. 
We are grateful to the referee for helpful comments. 
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Table 1. Spearman Rank-Order Correlation Coefficients 



Plot name 


Figure number 


Correlation Coefficient 


Two-sided significance 


L x vs L m 


El 


0.63 


5.0 x 10~ 5 


F x vs F m 




0.14 


0.40 


L x vs 8/5 


3 


0.38 


3 x 10- 3 


LlAGHz vs L IR 


4 


0.41 


0.02 


Fi.4GHz VS Fir 




-0.09 


0.61 


L Ha vs Lir 


5 


0.91 


3.6 x 10~ 12 


F Hq vs F m 




0.65 


1.1 x 10~ 4 


L(Ha) vs L(24//m) 




0.84 


2 x 10~ 15 


F(Ha) vs F(24/im) 


6 


0.67 


4 x 10~ 8 


M(H 2 ) vs Lir 


7 


0.95 


1.3 x 10~ 10 


F(CO) vs Fir 




0.81 


1.7 x 10~ 5 



Note. - Col 1. Correlation being tested. Col 2. Figure which plots the data. Col 
3. Spearman Rank-Order Correlation Coefficients. Col 4. Two-sided significance of the 
correlation's deviation from zero. The most signficant correlations are that between H Q and 
infrared luminosity and that between molecular gas mass and infrared luminosity. Most 
correlations are done both on fluxes and luminosities. 
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Table 2. Star Formation Rate(SFR) 



Cluster 




SFR 




/ n A A ~1 \ 

(10 44 erg s" 1 ) 


(M /yr) 


Zjzuoy 


R4 R8 


971 

Z 1 1 


A91 zLfi* 

r\Z 14fcU 


4^ 4fi 


1 Q9 
iyz 


/\1UOO 


44 fil 


1 88 
loo 


R0821 -4-07* 


8 47 

O . T: 1 


X! 


ja,iooz-|-ou 


99 R9 
zz.uz 


Q7 
y i 


781 

Zjoiyo 


i 3 7n 

id. 1 u 


oy 




1 1 Q9 

i i.yz 


^9 

OZ 


Anm i * 


7 Q7 


oo 


PKS0745-1 

-L 1VUU 1 TT'.y X 


3 8D 


1 7 2 


Al fifi4 


^ 91 


1 4 

14:. U 


Pfl^9-Ul Q 


9 4fl 


1 1 1 
11.1 




u.ou 


a n 


Rrms-i-nQ 




9 1 

Z. 1 


Dfi/i QQj_nc;* 
ruj^oy-ruo 


4 1 7 


1 8 7 
lo. ( 


A 9904 


3 9^ 


14 7 


A 9R97 


1 

i . oy 


7 ^ 


Am 1 ^ 


1 


9 
u.z 


781 Q7 


n 79 

U. 1 z 


o.o 


R9i 9Q_i_nn 
rtzizy-i-uu 


9 QQ 

z.yo 


1 Q /I 


A 1 9D/1 


1 7*? 

1. Id 


8 1 
o. 1 


A0646 


1.49 


7.1 


A2055 


1.46 


7.0 


A0291 


1.30 


6.3 


A1885 


1.04 


5.1 


A3112 


0.84 


4.2 


A2292 


0.80 


4.0 


A1930 


0.75 


3.8 


Z8276 


0.74 


3.7 


A4095 


0.29 


1.6 


A0085 


0.28 


1.6 
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Table 2 — Continued 



Cluster 


(10 44 erg s- 1 ) 


SFR 

(M /yr) 


A2052 


0.24 


1.4 


R0000+08 


0.20 


1.2 


NGC6338 


0.18 


1.0 


R0751+50 


0.10 


0.65 


A0262 


0.08 


0.54 



Note. Infrared luminosities are 



from iQuillen et al.l (120071 ) estimated 



from the 15/xm wavelength for BCGs 
that are detected at 70/^m or have color 
ratio F 8f , m /F 5 

.8/Ltm > 1.0 Or i*24/im/-^8/im ^* 

1.0. The star formation rate was esti- 
mated using Equation [TJ The top sec- 
tion contains four BCGs that are sus- 
pected to harbor dusty Type II AGNs. 
Z2089, A2146 and A1068 exhibit a red 
^4.5/Ltm/-p3.6/Lim color and all four exhibit 
high [OIII](5007)/H/3 flux ratios. Note 
that if there is an AGN present in one of 
these clusters, the SFR may be overes- 
timated. The second set is the remain- 
ing 10 BCGs with iWF 5 . 8Atm > 1.3. 
The third section is the set of 6 clus- 
ters with 1.0 < F 8 ^ m /F 5 , 8fim < 1.3. The 
fourth set is the remaining BCGs with 
IR excesses. Specifically they have ratios 

Fsfim/ F^^nm > 1-0, F2^ m / Fs.o^iia > 1-0, 

or a detected 70/xm flux. The BCGs 
marked with a * can be classified as 
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LIRGs since they have Ljr greater than 
lO n L . Objects with F 8Mm /F 5 . 8Mm < 1.0 
or F24, fim ./ Fsfaa < 1 and lacking a 70 /xm 
de tection are listed in Table 3 of paper 
I (j Guillen et al. 2007 ) with upper limits 
on L IR . For these objects L IR < 0.3 x 
10 44 erg s _1 and corresponding star for- 
mation rates are lower than < 1M Q yr _1 . 



Table 3. Comparison of Estimates of Star Formation Rate 



BCG 


OD07 


C99 


HM05 MQ93 


M95 


M04 


M05 


B03 M089 OD04 D07 




U.O 






0.02 


















12 








1 n a 

1U 4 


A1795 




2 


9 


12 








1.8 10 


A1835 




77-125 


123 








100 




55 

HydA 

55 






9.5 




1(b) 
23-35 (c) 




138 (FIR) 




A2052 


1.4 


0.96 












0.4-0.6 0.16 


A1068 


188 


30 








16-40 






A1664 


14 


23 








68 (IR) 






R1532 


97 


12 














A2146 


192 


5.6 














PKS0745 


17 




129 













Note. — Comparison of star formation rates from this paper (Col 2) with estimates from the literature (Col 3-12) Values 
arc in units of Mq yr _1 . References: OD07 = this paper, C99 = Crawford et al. (1999), HM05 = Hicks & Mushotzky (2005), 
M093 = McNamara & O'Conncll (1993), M95 = McNamara (1995), M04 = McNamara et al. (2004), M05 = McNamara ct al. 
(2005), B03 = Blanton ct al. (2003), M089 = McNamara & O'Conncll (1989), OD04 = O'Dea ct al. (2004), D07 = Donahue 
ct al. (2007). For Hydra A, (b) and (c) refer to short duration burst and continuous star formation models, respectively. 
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Table 4. Mass Deposition Rates 



uluster 


M s 


Mj 




(M /yr) 


(M /yr) 


AlUDo 


qn + 20 

<3U_ 10 


44U_ 10 


R1532+30 


4ooii88 


i9ooiq 


PKS0745-1 


200^8 


10801™ 


A1664 


60+1° 


3301^8 


R0338+09 




27018 


A2204 


' u -40 


860188 


A115 


°-6 


190118 


R2129+00 


6lf 


3801^8 


A1204 




620118 


A3112 


101* 


220118 


A4059 


5+ 2 
d -i 


105+3 


A0085 


6±S 


210II8 


A2052 




7211 


A0262 


1-818.1 


10+1 



Note. - Mass deposition rates 
calculated within r coo i using Chandra 
and XMM data. M s is a measure of 
the mass deposition rate consistent 
with the X-ray spectra and Mj is a 
measure of the mass deposition rate 
if heating is absent. Sources are listed 
in order of decreasing SFR. 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Redshift 



Fig. 1. — Flux ratio F4.5/F3.6 vs. redshift. The labeled objects show evidence for the presence 
of an optically luminous type 2 AGN. 
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1e+45 



CO 

x 



1e+44 



1e+43 



<- 
% 4 , 



^2089 

(A2146 
1 A1068 



^ 



aA2292 
H M 



N4104 



■R0821+07 
IA11 



F 8 /F 5 o > 1 .3 
1.0<F 8 /F 5 r<1.3 

«8 l %8 <1 ;f24 /F 8 >1 

Additional F 70 



▲ 
♦ 



1e+43 



1e+44 1e+45 
Lir (erg/s) 



1e+46 



Fig. 2. — X-ray luminosity vs. estimated infrared luminosity. X-ray and infrared luminosities 
are listed in Tables 1 and 3 of Paper I. Solid squares are our reddest group with 8 to 5.8/im 
flux ratio greater than 1.3. The intermediate group with flux ratios between 1.0 and 1.3 are 
plotted as stars. Solid triangles have 8 to 5.8 /zm flux ratio less than 1 but 24 to 8 fim flux 
ratio above 1. Solid diamonds are galaxies with both 24 to 8 and 8 to 5.8 /im flux ratios 
less than 1 but have been detected at 70/xm. Flux ratios are computed using photometry 
listed in Table 2 in Paper I. Upper limits on the infrared luminosity are shown by arrows. 
We find a weak correlation between X-ray and infrared luminosity. The dashed line is the 
kinetic energy injection rate predicted from a star forming population due to supernovae as a 
function of the infrared luminosity. We confirm that kinetic energy from supernovae cannot 
account for the X-ray radiative energy losses in most cooling flows. 
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CD 



1e+45 



>< 1 e+44 



1e+43 



□ R1 532+30 




0.7 



6 DA1068 



□ 



□ 



^A1664 
^□R0338+09 



□A11 



R0821+0^ 



□ N4104 



1 2 
F(8|im)/F(5.8|im) 



Fig. 3. — X-ray luminosity vs. 8 to 5.8/xm flux ratio (data taken from Tables 1 and 2 in 
Paper I). Most red objects with F 8fl / F 5 . 8fim > 1 have X-ray luminosity L x > 10 44 erg/s. 
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1e+28 
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R0821+07 
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8' 1 5.8 



1.3 



IN4104 1 -0 < F 8 /F 5^ < 1 
F Q /Fco<1, f 0/1 /f 



ional F 



.' '24" 8 



70 



3 * 
>1 A 

♦ 



1e+43 



1e+44 



1e+45 



1e+46 



L|r (erg/s) 



Fig. 4. — Radio luminosity at 1.4 GHz (listed in Table 1 of Paper I) vs. estimated infrared 
luminosity (Table [2j The radio- infrared correlation (equation [2]) for star forming objects is 
shown as a line on the lower right. The radio fluxes are much higher than expected from the 
radio-IR correlation appropriate for star forming late type galaxies. This is not unexpected 
since many of these objects contain radio cores and in some cases even double lobed jets. We 
find a weak correlation between the radio luminosity at 1.4 GHz and the infrared luminosity. 
The point types are as given in Figure [21 
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Fig. 5. — Observed Ha Luminosity (listed in Table 1 of Paper I) vs Infrared Luminosity 
estimated from the 8/xm and 24.4/im fluxes. The data for 2 BCGs from lEgami et al.l (j2006al ) 
are shown as open circles. The Kennicutt relation inferred from observations of star forming 
galaxies relating Ha luminosity to star formation rate is plotted as a line. We have divided 
the line by a factor of 2.8 to back-out the reddening correction since our Ha luminosities 
are uncorrected for reddening. The Ha fluxes are consistent with the estimated infrared 
luminosities and star formation. As is true for Figure [2] the point types depend on the 8 to 
5.8/iin color. We suspect that some of the Ha luminosities are lower than expected because 
the apertures used to measure them were smaller than those used to measure the infrared 
fluxes. 
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Fig. 6. — Observed Ha flux vs 24 /iin flux listed in Paper I. Arrows denote upper limits. 
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Fi g. 7. — Correlation o f H 2 mass from CO lumin osity and Ltr. A1835 a nd Z3146, discussed 



by lEgami et al.l (j2006al ). and A2597 discussed by lDonahue et al.l ( I2007bl ) are shown as filled 
circles. We find a strong correlation and as such consider the relation between CO and star 
formation. Upper limits are shown as arrows. The 4 objects thought to host AGNs are 
labeled. Two of these have IR luminosity higher than expected from their molecular gas 



mass. 
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Fig. 8. — Co rrelation of mass f rom CO luminosity and SFR. 2 BCGs, A1835 and Z3146, 
discussed by lEgami et al.l (I2006ar) . are shown as filled circles. The dotted lines are taken 
from Eq. 4 by iKennicutt Jl998l ) using; different values for the diameter of the star forming 
region. See the legend at the upper left. The diameter of the star forming region tends to 
be larger for more luminous objects. Solid lines are computed assuming gas depletion time 
scales of 10 8,5 , 10 9 and 10 9 5 years. See the legend at the upper left. 
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Fig. 9. — IR luminosity vs. X-ray derived cooling times at a radius of 30 kpc. BCGs have 
higher IR luminosity in clusters with shorter cooling times. 
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Fig. 10. — X-ray derived mass deposition rate upper limits against estimated star formation 
rates. The closed circles correspond to maximum mass deposition rates, Mj if heating is 
absent, and the open circles refer to M s , the mass deposition rate consistent with the X-ray 
spectra. The dashed line is for equal star formation and mass deposition rates. 
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Fig. 11. — The X-ray luminosity emitted with the cooling radius (a fiducial radius where 
the radiative cooling time is 7.7 Gyr, corresponding to redshift one) is plotted in solid circles 
against the Spitzer infrared luminosity. The expected (missing) luminosity emitted below 
1 keV by a continuous cooling flow operating from the cluster virial temperature to zero K 
is shown by the open circles. If mixing with dusty cold gas causes the rapid non-radiative 
cooling of the intracluster gas below 1 keV then this luminsoity could emerge in the mid- 
infrared. 



